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Abstract 

In various background independent approaches, quantum gravity is defined in terms of a field 
propagation kernel: a sum over paths interpreted as a transition amplitude between 3-geometries, 
expected to project quantum states of the geometry on the solutions of the Wheeler-DeWitt 
equation. We study the relation between this formalism and conventional quantum field theory 
methods. We consider the propagation kernel of 4d Lorentzian general relativity in the temporal 
gauge, defined by a conventional formal Feynman path integral, gauge fixed a la Fadeev-Popov. If 
space is compact, this turns out to depend only on the initial and final 3-geometries, while in the 
asymptotically flat case it depends also on the asymptotic proper time. We compute the explicit 
form of this kernel at first order around flat space, and show that it projects on the solutions of all 
quantum constraints, including the Wheeler-DeWitt equation, and yields the correct vacuum and 
n-graviton states. We also illustrate how the Newtonian interaction is coded into the propagation 
kernel, a key open issue in the spinfoam approach. 



1 Introduction 

The tentative quantum theories of gravity that are currently better developed, such as for instance 
strings and loops, are very different from one another in their assumptions and in the formalism 
utilized. In particular, the relation between the background-independent methods used in canonical 
quantum gravity and in the spinfoam formalism, and conventional perturbative quantum field theory 
(QFT), is far from transparent Besides clouding the communication between research communities, 
these differences hinder the clarification of a number of technical and conceptual problems. For 
instance, there are open questions concerning the physical interpretation of the spinfoam formalism 

2 , in particular concerning its low-energy limit, the comparison with quantities computed in the 
standard QFT perturbative expansion, and the derivation of the Newtonian interaction. 

Here we contribute to the effort of bridging between different languages, by studying the field prop- 
agator, or Feynman propagation kernel, or Schrodinger functional, in 4d Lorentzian general relativity 
(GR). The field propagator (and its relativistic extension PP) is not often utilized in conventional 
QFT (but see [H] and 0j), but it is a central object in background-independent quantum gravity 
El El Ej ■ Here we analyze this object starting from a conventional path integral quantization of GR 
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and -neglecting its non renormalizability- in a perturbative expansion of this integral. On the one 
hand, this provides a clean interpretation of a basic non-perturbative tool in terms of conventional 
and well understood quantum field theoretical quantities. On the other hand, this provides the precise 
expression of the low-energy limit of the field propagator, to which the non-perturbative one must 
be compared. In other words, we study the explicit relation between the 3-geometry to 3-geometry 
transition amplitude and the graviton-state language. 

The key to bridge between a conventional path integral formulation and the non-perturbative 
framework is the use of the temporal gauge, with a careful implementation of the Faddeev-Popov 
(FP) gauge-fixing. Conventional perturbation theory is usually studied in covariant gauges such as 
the Lorentz gauge in Yang-Mills (YM) theory or the harmonic gauge in GR ^UJi but the temporal 
gauge is naturally closer to the Hamiltonian formalism. The propagation kernel of YM theory in the 
temporal gauge is well understood It can be computed as a Feynman path integral for finite 

time, with fixed initial and final field configurations. In the temporal gauge there are no dynamical 
ghosts, and the kernel has a clear physical interpretation: it gives the matrix element of the evolution 
operator between eigenstates of the YM connection. Namely, the transition amplitude between the 
states defined on the boundaries. But it is also a projector on the physical states, which satisfy the 
YM quantum constraints. At the zero'th order in perturbation theory, the kernel can be explicitly 
computed with a gaussian integration in the Euclidean regime, and it nicely codes the form of the 
perturbative vacuum as well as all the n-particle states Furthermore, it is straightforward 

to express the n-point functions in terms of it. In spite of key differences, the structure of GR is similar 
to a YM theory in many respects. We can apply to GR the techniques used for YM theory, and, in 
particular, study the propagation kernel of quantum GR in the temporal gauge. This is what we do 
in this paper. Notice that the use of non-covariant gauges has been considered in the literature (see 
for instance [T3|). 

First, we consider the formal path integral that defines the propagation kernel in the temporal 
gauge, and study its properties. This object has been considered in the literature (see |H1 El and 
references therein), though with different techniques. The expression is formal because the measure 
is unknown, and the usual perturbative definition is not viable because the perturbative expansion 
on a background is non-renormalizable. We do not consider background-independent definitions of 
the path integral, such as the spinfoam one, because our interest here is the interpretation and the 
low energy limit of the kernel, not its ultraviolet divergences. We consider the two possibilities of 
compact and asymptotically-flat space JS]. We carefully discuss the FP gauge fixing and in this 
context we identify the integration implementing the quantum constraints as the one on the FP gauge 
parameters. We discuss how the propagation kernel turns out to be independent from the coordinate 
time in both cases -a feature that drastically distinguishes GR from YM theory-, but to depend on 
the asymptotic proper-time in the second case. We use the conventional metric formalism for GR. 
The expression we obtain is formal, and we do not discuss topological and ultraviolet aspects. 

Second, we consider the zero'th order term in the perturbative expansion of the integral on a flat 
background, and we compute the propagation kernel explicitly. We show that it projects on the solu- 
tions of all the constraints, and that it correctly codes the perturbative GR vacuum state |16M17| and 
the n-particle states. This provides an explicit bridge between 3-geometry transition amplitudes and 
perturbative graviton states. The non-perturbative boundary amplitudes of the spinfoam formalism 
must reduce to this expression for boundary metrics close to flat space. 1 

Finally we couple an external matter source to the theory and derive the expression for the energy 
of the field in the presence of matter. This expression codes the Newtonian interaction. Our hope is 
that this could open the way for extracting the Newtonian interaction from the spinfoam amplitudes, 



1 For recent applications of this idea, see |18l . 



2 



hence providing a key missing check of their physical viability. 

We fix the speed of light and the Newton and Plank constants by c = WirG = K = 1. Greek indices 
range from to 3; latin indices from 1 to 3. We use coordinates x^ — (t,x). A fJi (x) is the YM field, 
Ai(x) the vector potential; g^ v {x) is the gravitational field and the spacetime metric; while gij{x) is 
the metric of a spacelike surface. (Indices are intended in Geroch's abstract index notation: A4 means 
A, and so on.) 



1.1 YM propagation kernel in temporal gauge 

Before turning to GR, we recall the properties of the field propagation kernel in YM theory. As we shall 
see, the gravitational case will present substantial analogies, as well as key differences. Formally, the 
propagation kernel is given by the functional integral over the configurations of the YM field defined 
on the spacetime region bounded by the initial and final surfaces t — and t = T, and restricted to 
given initial and final configurations A^x, 0) = A' (x) and A^{x, T) — A''(x) 



W[A'^A>> 



T] = f, VA^ e iS ^\ (1) 



where 

fT 



(i 



S[A^ = dt d 3 x C YM (2) 



and £ym is the YM lagrangian. The integral JQ) contains an infinity due to the integration over the 
group of the gauge trasformations A — > A A , where A(x, t) is the gauge parameter. We fix this by 
gauge fixing Q a la FP in the temporal gauge Aq = That is, we insert in Q the identity 

1 = AMA»)[VA6(A£). (3) 

The FP determinant App is actually a constant, because the YM gauge transformations do not mix 
the different components of the 4-vector A^, hence App can only depend on Aq] but this, in turn, is 
fixed to be zero by the ^-function appearing in the integral. The propagation kernel between boundary 
values of the 3d YM connection Ai is thus given by 

W[A!i,A!l,T\ - App J aH VA^j VAS(A^) e tS ^l (4) 

A 1 . 

Changing the order of the two integrations, changing variables A^ — > A A and integrating over Aq, we 
obtain 

WW, Al T] = App f VA f VA t e ^A =o]_ (5) 

The only components of the integration variable A entering the integrand are its values at t = and 
t = T. We can therefore drop the bulk integration on A(x, t), > t > T, discarding a trivial infinity. 
Furthermore, since VAi and S[Ai, Aq = 0] are invariant under time independent gauge trasformations, 
the T>Ai one of the two remaining integrals on A(x, 0) and A(x, T) is redundand. Dropping the second, 
we have therefore 

W[A' i: A", T}= Jv\ WIA'S, A'l, T] (6) 

where 

W[A^,AH,T\ = App f VAi e iS ^ Ai ' Ao=0 \ (7) 
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and X(x) is the gauge parameter of the residual time-independent gauge trasformations Ai — ► A\ of 
the Aq — gauge. The propagator W[.A£, A'( , T] is invariant under simultaneous gauge trasformations 
on the two boundaries. The integral over A in © makes W\A^A",T\ invariant under independent 
gauge trasformations on the two boundaries. 

The field propagator is the matrix element of the evolution operator between eigenstates of the 
field operator A.- L 

W[A' i ,A'l,T} = (A!>\e- iHT \A! i ), (8) 

where H is the hamiltonian; up to the difficulties in defining fixed-time operators in an interacting 
QFT (see [3]), it can be interpreted as the Feynman probability amplitude of having the field config- 
uration A" at time T, given the field configuration A i at time 0. Equivalently, it time-propagates the 
quantum state in the Schrodinger functional representation of the quantum field theory 

* t+T [M = J VA\ WlAuA'^T] %[A'A. (9) 

However, notice that, because of its gauge-invariance, any state obtained by propagating with the 
kernel is invariant under gauge trasformations of Aj. Therefore the kernel is also a projector on the 
gauge-invariant, or "physical" , states, which satisfy ^[Ai] — ^f[A^]. That is, the states that satisfy the 
Gauss-law quantum constraint. Thus, the T>Ai integral (JJJ takes care of the gauge- variant dynamics, 
and the T>\ integral © imposes gauge invariance. 

If ^ n (Ai) is a basis of gauge-invariant physical states that diagonalize the energy, (JSJ implies 

W[A h A' il T]=Y,e- iE - T %JM^[A i ]. (10) 

n 

In particular, once subtracted the zero-point energy, we can read out the form of the vacuum state 
from the propagator 

/+oo 
dT W[Ai,A'i =0,T} = const * [^i]- (H) 
-oo 

Thus, the temporal-gauge propagation kernel nicely bridges between the functional integral formalism 
and the hamiltonian one. 

The propagation kernel W\A! i ,A", T] can be computed explicitely order by order in perturbation 
theory, by a gaussian integration in the Euclidean regime. For this to be well defined, we need to 
assume appropriate boundary conditions for the field at spacial infinity. In particular, we assume that 
the field vanishes at spacial infinity, and therefore so has to do the gauge parameter. In the lowest 
order (or exactly in the Maxwell case) the propagator kernel turns out to be [TT] 

d 3 p (\A' T \ 2 + \A" T \ 2 )co S pT-2A' T -A" T 

WiA^A^T] = Af(T)e/J (27r)3P sinpT . (12) 



Here p — \p\ — y/ Pip 1 , the Fourier trasform of the potential is Aj(p) — J d 3 x e l ^' x Aj(x) and its 
transverse component is defined as Af(p) = (p)Aj(p), where 

D l3 =5l-?f. (13) 

We can read out from this kernel all the free n-particle states with momenta p a , energy E a = p® a and 
polarizations e Q , a = 0, ...,n, from the expression 

W[A i ,A> i ,T] = ±± E / ^ • - TO e-^ E « T i^ZZZM] ^,..., Pn M (14) 

' n=0 ei...e„ J ^ ' ^ ' 

Let us now turn to general relativity. 
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2 Propagation kernel in general relativity 



We foliate spacetime with a family of 3d surfaces Ej, with t € R, and and focus on the region t G [0, T]. 
We fix initial and final positive definite metrics g'^ on So and g'/j on Et as boundary data, and we 
want to describe the quantum dynamics of the gravitational field in terms of the functional integral 

W[g' ij ,g^T} = J S{9 »» ] (15) 

over the 4d spacetime metrics inducing the given 3d metrics on the two boundaries. The action of the 
gravitational field is 

S[g„ v ] = [ dt [ d 3 x y/^jg^R^ + [ d 3 x K. = [ dt [ d 3 x C. (16) 
Jo JT, t JsouSt Jo J'St 

Here R^ is the Ricci tensor and g the determinant of the metric. K. is the extrinsic curvature of the 
boundary surfaces. The presence of the boundary term, sometimes called the Gibbons-Hawking term, 
is needed in order to have only first-order time derivatives in the action C, so that the convolution 
property of the propagation kernel is guaranteed |19) . 

There are various sources of infinities in (|15|l . First, there are ultraviolet divergences. As dis- 
cussed in the introduction, we disregard them here, under the assumption that an appropriate non- 
perturbative definition of the integral, such as in the spinfoam formalism, could take care of them. 

Second, we must be sure that a sufficient number of boundary conditions are fixed. In general, we 
may reasonably demand that a classical solution is uniquely selected by these conditions by minimizing 
the action. We focus on two cases: (i) The compact case, in which the Riemaniann manifolds (Eo, g[j) 
and (YiT,g'ij) have finite volume and no boundary. In this case, initial and final 3-metrics may 
determine a unique classical solution. The thin sandwich conjecture [HI I2()|. indeed, states that, at 
least for small times, generically there is only one spacetime metric between two given spacelike 
metrics. Exceptions are known, generally characterized by pathologies such as singularities [22- In 
the following we restrict to the cases where the conjecture holds. Amongst the exceptions is the 
remarkable case of flat space, as emphasized in [32]. (ii) The asymptotically flat case, in which E is 
homeomorphic to R 3 and g'^ and g'/j appropriately converge to &y at infinity. In this case, for 115(1 
to be well defined, we demand g^ v to converge asymptotically to the Minkowski metric rj^ for all t. 
The consequences for the uniqueness of the solution will be discussed in details below. 2 

Third, the invariance under diffcomorphisms of GR makes the integral i|15|) infinite. This situation 
is analogous to the YM case, and can be cured with a gauge-fixing, as we do below. 

2.1 Gauge— fixing 

General relativity is invariant under diffeomorphisms, namely under the pull back g — > g^ = 

of the gravitational field by a map £: M — > M from the spacetime to itself. These are the gauge 

transformations of GR. 3 Explicitly, g^ is defined by 

9 ^ {x) ^^y^x^ 9 ^ {x)) (i7) 



2 It has been recently suggested that the most interesting case for nonperturbative quantum gravity is when St has 
boundaries 1231 . For a discussion of the classical solutions in this case, see 1241 . 

3 They should not be confused with the freedom of chosing coordinates on M: once coordinates are fixed, the integral 
1151 still gets contributions from distinct fields g and g^. 
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where £ : x^ i— > ^(x). (Or g^(x) = fjjrf^ g pa {(j)(x)) where (/) = a form we'll use later on.) For 
the action (|16J) to be invariant, £ must not change the boundaries of the spacetime region considered, 
that is 

£>,£) = 0, e{T,x)=T. (18) 

Because of this gauge invariance, the integral (|15JI has an infinite contribution from the integration 
over the gauge group. We take care of this as we did for YM, by introducing a (non-covariant) gauge- 
fixing. The GR analogue of the temporal gauge is known as "gaussian normal coordinates", or the 
"Lapse=l, Shift=0", or "proper-time" gauge 

3oo = -1, goi = 0. (19) 

As for YM, this is not a complete gauge-fixing, but we expect that additional gauge-fixing is not 
required in the path integral. In the linearized case we shall explicitly see that the remaining part of 
the gauge is taken care by the integration over the gauge parameters. Thus, we gauge-fix the path 
integral by inserting in (|15fl the FP identity 

1 = Appfo^] J Vi %« + 1) 5(gi). (20) 

T>t; is a formal measure over the group of the 4d diffeomorphisms. We can see here a first difference 
with YM theory: the GR gauge transformations mix the different components of g pu , hence the FP 
factor App is not a constant anymore. Since the four go^ are fixed by the ^-functions, it will depend 
only on the spacial components of the metric. 

The integration Il20|) is over all £^(x, t) with t G [0,T], including the boundaries, and it is not 
restricted by Ijl8|l . therefore it includes £ that change the action (|16|l . To understand this delicate 
point, observe that the FP integral must include sufficient gauge trasformations for transforming any 
field to one satisfying (|19fl .We cannot fix goo = — 1, and also the coordinate time between initial and 
final surface. This would amount to discard all four-metrics yielding a proper time between the two 
surfaces different from T; but these field configurations do contribute to l|15|) and cannot be discarded. 
In other words, we cannot gauge transform all fields contributing to the integral to the gauge l|19|) 
without changing the action. However, nothing prevents us from transforming them to fields satisfying 
i|19fl and changing the action when needed. (See also [S] on this.) 

Introducing (|20|l in (|15J) we obtain 

W\g , ij ,g? j ,T] = J g „ Vg^ArpfaA J V£ exp J d 3 x £ dt £m|%cSo + 1) %tk)- ( 21 ) 

We now evaluate the integrals over go^. As we did in the YM case, we exchange the order of the 
integrations and change variables g^ u — > g^, together with a change of coordinates x M — * £ M (x) in the 
action. As a consequence, the boundary data will now depend on £. We perform the integrals over 
goo and g i, obtaining 



W\g' ip g'l j ,T} = Jvs, U K.^.V (22) 



where 



Wlg'i^g'^t, T] = I t Vg^ A PP [p tf ] exp | i [ d 3 x [ dt C[g VJ ,goo = l,5oi = 0, ] 1 . (23) 



Let us analyze these expressions. The T>gij integral in l|23|) is over the fields in the spacetime region 
bounded by the surfaces Eq and Ej., defined respectively by t — £°(x, 0) and t — £°(x, T), having 



G 



boundary value g"£ on Eg and g[£ on Ey. Notice that gL^ depends precisely on g[- and £, since, 
using (fli^l , the transformation of gij reads 



s a? a? oe deeded? 



that is, the value of g[£ on Eg is determined by on Eo and by the map £ : Eo — + Eg. 

Equation (|23[1 depends only on the boundary diffeomorphisms Cku(^) = £ m (£j0) an< i CfinO^) = 
^(x,T). The integral over the bulk diffeomorphisms depends in general on the boundary surfaces. 4 
The result of this integration is an appropriate functional such that the convolution property of the 
kernel will be satisfied. However, since we are interested in explicitly computing the linear approx- 
imation, this term is not relevant and it will be discarded in the following. The boundary gauge 
transformations play a more subtle role than in YM theory. Indeed, notice that we cannot write 
Wlg'ij, g"j,£,T] in the form W\g'^,g'^ ,T] as we did for YM, because £ affects also the action. Let us 
distinguish the temporal boundary diffeomorphisms and £g n from the spacial ones, £? ni and £ fin . 
Truly, (12.SII is invariant under a spacial diffeomorphism which acts identically on both boundaries. This 
fact is analougous to the YM case, and allows us to drop the dependence on £* of one of the boundaries, 
say at t = 0. On the other hand, this property is not true for temporal diffeomorphisms as emphasized 
in 0. Thus we can write, introducing the shorthand notation = C[gij,goa = ~l,3oi = 0], 

H '/;,.//;;. / - J nL j g ,, r v 9lJ a kp ^] exp u J d 3 x dtc[ g A. (25) 

q'. . £® ^ ^ 

This is the exact expression for the propagation kernel of GR. We expect the integration over the 
gauge parameters to implement the constraints of the theory. In section 3 below, we show explicitly 
how this happens in the linearized case. Before that, however, we discuss a feature of the kernel which 
is characteristic of GR. 

2.2 The disappearance of time 

Let us first focus on the case (i) in which (E, g^j) is a compact space with finite volume, and consider 
the 2?£g n integration in This is an integration over all possible coordinate positions of the final 

surface. This ensemble is the same whatever is T, hence the right hand side of l|25|l does not depend 
on T. Therefore in this case the propagation kernel of GR is independent from T, and we can simply 
write 

H '/;,•//" = J 2>&, 2>& V&ifw V m AFp[ff -]exp|i J d 3 xj^dtC[ 9lJ ]^. (26) 

This is what we mean by disappearance of time. Notice however that a proper time is (generically) 
determined by the initial and final 3-metrics themselves. In fact, the integral (|26|l is likely to be picked 
on the classical solution <? M „ bounded by g'^ and g'[y But recall that in classical GR the initial and 
final 3-metrics are expected to generically determine a classical time lapse between them. To see how 
this can happen, consider the theory in the partial gauge-fixing <?oi — 0- The six evolution equations 
for gij (x, i) are second order equations that we expect to generically admit a solution for given initial 



We thank our referee for pointing this out. 
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and final data at t — and t = 1. The time-time component of the Einstein equations -the scalar 
constraint- can then be written in the form 



//V - r/ J 'i 



) 



dgjj dg k i 
dt dt 



goo det gij R[g i:j } = 0, 



(27) 



where R is the Ricci scalar of the 3-metric gtj . Once a solution g^ [x, t) is given, this equation can be 
immediately solved algebraically for goo (x, t) and therefore it determine the physical proper time 



between initial and final surface, along the x = const lines, which in this gauge are geodesies normal 
to the initial surface. Hence in general this proper time is determined by the initial and final 3 
geometries. On the other hand, notice that as an equation for goo, (|27() becomes indeterminate when 
R[gij] — 0, and in particular on flat space. 

The disappearance of the time coordinate in (|26J) . and in general in the transition amplitudes of 
quantum gravity, and its physical interpretation, have been amply discussed in the literature. Its 
physical meaning is that GR does not describe physical evolution with respect to a time variable 
representing an external clock, but rather the relative evolution of an ensemble of partial observables. 
See for instance for a detailed discussion, and the references therein. 

The invariance of GR under coordinate time reparametrization implies that the hamiltonian H in 
JSJ vanishes, and therefore we expect to have 



instead of (|10(l . where again ^[g^] is a complete basis of physical states, satisfying all the Dirac 
constraints of GR, including, in particular, the Wheeler-DeWitt equation, which codes the quantum 
dynamics of the theory in the hamiltonian framework. Equation (|29(l indicates that W[g' i: j, g'/j] is the 
kernel of a projector on the physical states of the theory j2HJ- Roughly speaking, the 2?£ l integration 
implements the invariance under spacial diffcomorphisms, while the 2?£° integrations project on the 
solutions of the Wheeler-DeWitt equation. We will see that this is indeed the case in the linear 
approximation. 

An explicit perturbative computation of the propagation kernel in the compact case would be very 
interesting. The problem of expanding around a flat background in the compact case is that this is 
precisely one of the degenerate cases where the thin-sandwhich conjecture fails; essentially because, 
as we have seen, (I27|l fails to determine goo in this case. However, this can probably be simply 
circumvented, for instance by adding a small cosmological constant. We leave this issue for further 
developements, and we turn, instead, to the asymptotically flat case, where the explicit computation 
of the propagation kernel can be performed in a more straightforward fashion. 

Thus, consider the case (ii) in which (E,^) is asymptically flat. In this case, again l|27(l does 
not determine T. However, the argument above for the disappearance of time fails, for the following 
reason. As we have already mentioned, for the path integral to be well defined we must require flat 
(i.e. g^ = 5ij) boundary conditions at spacial infinity for gij(x,t). In order to be well defined on 
this space of fields, the gauge trasformations must vanish at infinity accordingly, as they are required 
to do in the YM case. Therefore £°(x,T) must converge to £°(x,T) = T for large x, and the right 
hand side of Ij25(l is not independent from T. In fact, what is physically relevant is obviously not the 
coordinate time, but rather the proper time separation between initial and final surfaces at infinity. 
That is, in the asymptotically flat case, the GR field propagator depends on the asymptotic proper 



T(x) = / y/g m (x, t) dt 



Jo 




(28) 




(29) 



n 
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time at infinity. 5 



3 Linearized Theory 

We now focus on the asymptotic case, write 

9nu{x) = t]^ + h^(x) (30) 

and consider the field h lJLll {x) as a perturbation. Since we restrict to small fields, we restrict, 
accordingly, to small diffeomorphisms, that preserve the form l|3U|) . The gauge condition (|19fl implies 
hopi(x) = 0. In this gauge, the action i fTfijl reads 

■T 

d 3 x I dt £ = I d 3 x I dt { didjh 13 - V 2 h 



o J Jo ^ 

+~ d a h\d a h\ - d a h ij d a h ij + i) k >>";),>,,,, + d k h ij dih jk - d l h ll d h ~ if !,, ,<)■ I, j . (31) 

Indices are now raised and lowered with the Minkowski metric and the 3d Euclidean metric. It is 
convenient to separate the four Poincare-invariant components of the linearized field hif. the spin- 
zero trace of the transverse components ho, the spin-one and spin-zero longitudinal components 
hi = hj + ^j-ft. L , and the spin-two traceless and transverse components hJJ . In Fourier space, this 



decomposition is 

where Dij is the projector on the transverse modes given in (|13|) and 



hi, = -2/iqA, + ^- + hj- + ^ L ^r + ^ T (32) 
P p p 



p i hg = 0, ^^/ = 0, ^ = 0. (33) 

See Appendix C for more details. 

The action is invariant under infinitesimal spatial diffeomorphisms, x % — > C( x ) — % l + £l ( ;E )i which 
induce the transformation hij(x) — ► h e ij{x) — hij(x) + di£j(x) + djti{x) on the linearized field. Or 

hi{x) -» h e i(x) = hi(x) + £i(x). (34) 

Note, from this last expression, that £i(x) is of the same order of hij(x). 

As mentioned, the role of the temporal infinitesimal diffcomorphism is subtle: the lagrangian is 
insensitive to them, but they affect the extrema in the t integral in the action. 6 The term didjh 1 ^ — 
~S/ 2 h = — 2V 2 /io appearing in the lagrangian is a boundary term, and does not contribute to the 
equations of motion, but it does contribute to the integral. It is the linear approximation to the 
left hand side of the time-time component of the Einstein equations (|27|) , namely of the GR scalar 
constraint. 



5 A suggestive way of understanding the presence of asymptotic time in the asymptoticaly flat case, based on the 
boundary interpretation of the time evolution developed in HI 1231 is the following. According to HI 1231 . the observable 
proper time in GR can be determined by the boundary value of the propagation kernel of a finite spacetime region, 
namely a 4d ball. In the limit in which the spacial dimensions of the region go to infinity, the boundary proper time 
converges to the asymptotic proper time. 

6 Condition 1181 for the invariance of the action now reads e°(x, 0) = 0, e°(x,T) = 0. 
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The boundary data hij(0, x) — h'^ix) and hij(T, x) = h'^(x), vanish at spacial infinity. Hence l|25|) 
becomes 

W\h' 3 ,h>> 3 ,T] = fvei n Vel n V^ ni Vh* A FP [^] e ^ 3l C^ if *^, (35) 



The 2?£g n integration implements the invariance (|34|l separately on initial and final data. It thus makes 
W^j, h'(j, T] independent from the spin-one and longitudinal spin-zero longitudinal components h[ 
and h'l . How about the integrations over the e°'s? Since the background metric is static, e° enters 
only the action's boundaries. 7 Taylor expanding the action we have 



T+e°(x,T) 

j- I dt £[hij, h Q)1 = 0] 

0+e°{x,0) 
T 



i f d3x f a dt \ d a hd a h - a h lJ a h VJ + d k h lj djh ik + k h lj l h ]k - 0,1,' '0,1, - O^h^h 



+ d 3 x I e°(T,x) {0,0,}!, , - V 2 /i) 



t=T 



e°(0,x) (didjhij ~V 2 h) 



t=o 



(36) 



Thus, the integration over the e°'s gives two <5-functions of the linear term OiOjh 1 ^ — V/i on the 
boundary data 



W[h' ip h'!pT] = Sfoo^i'-V 2 ^") Sfah'y -V 2 //) fve 



J hlltl Vh i:i A FP (/iy ) exp I ~ / d 3 x J dt(0 a h0 a 

h'. . I 



h-0 a h i:> a h lit 



+O k h lj d 3 h ik + (/'//'•' (7, h ,i, - 0,1,' * 0,1, - O'h.jO'h 



(37) 



The argument of the first <5-function does not depend on e 1 , since, from the tranformation properties of 
hij we have did 3 h\" — \7 2 h e " = didjh'-j — S7 2 h". Recalling also that the constant Fourier component 
of the boundary data vanishes because of the conditions at infinity, and absorbing a constant in 
the normalization factor, we can write the 5-functions simply as 5 (/z )<5(/iq ) . These two 5-functions 
impose the linearized scalar constraint h$ = on the boundary data, and project on the solution 
of the quantum scalar constraint, namely the Wheeler-DeWitt equation. As emphasized by Kuchaf 
in [16| . indeed, in the linear approximation the scalar constraint is not anymore a relation between 
momenta and configuration variables, but rather a condition on the configuration space: the solution 
of the linearized Wheeler-DeWitt equation are the wave functions with support on ho — 0. Using this 
argument, in Hartle imposes these two 5-functions on the ground state functional for Euclidean 
quantum gravity by hand. Here, instead, we obtain them as a result of the FP integration calculations. 

Notice the two different ways in which the linearized temporal and spacial diffeomorphisms elim- 
inate, respectively, the spin-zero ho and spin-one and spin-zero longitudinal components hi of the 
lineraized field hij. The physical states are independent from hi, while they are concentrated on 
ho = 0. Finally, we are left only with the physical spin-two field hj^. 



7 On a generic background g9., the transformation under diffeomorphisms of the perturbation is, at first order, 
hij(x) — > h t ij(x) = hij(x) + di€j(x) + djei(x) + t^d^g^^x). Therefore, e° would enter the boundary data. 
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The FP factor App is cubic in the fields jOJ . Since in the following we are interested in computing 
the gaussian approximation, we neglect it. The gaussian integral over hij is straightforward; we give 
the details of the integration in the appendix. Writing 

H^(p) = W^Wh^W + W^h^j®, (38) 
the result of the integration is 

W[h , ij ,h'lvT]=M{T)8{h' )5{h'l)e li ^ P , (39) 

where the function Af(T) is a normalization factor. This is the field propagation kernel of linearized 
GR. 



3.1 Ground-state and graviton states 

We are now ready to read the vacuum state and the ro-graviton states from the propagation kernel 
H39fl . To do so, we expand (|39|l in power series of e~ lpT . To first order, we obtain 



W[h'h'li,T] = N{T) 5{h' Q ) 6(K) e-* J H (p) 



1 



1 



d 3 p 



P ff TT (p) e 



-ipT 



(40) 



Using (| 1 1|) . the (non-normalized) vacuum state can be read from the zero'th order of (|4t)(l : we have 
^\h%\MKj] = S(h' ) 5{K) exp { - \ J ^ P H^(p)}, 



and therefore 



Vofoj] = S(ho) exp{-J| ^ph^ l3 (p)h TT v(-P)} 



(2tt) 3 

This is in agreement with the literature |16l I17| . 

The graviton states can be obtained from the analog of (|14|l . namely 



(41) 



W \h- h> Tl-lV V f-^- ^e-^ 
W[hij,h ir l\ n ,2^ 2^ J (2^)3 ■ ' • ( 27 r)3 6 



E m T 



i>Piei,... 



n—0 ei...e 



This expression can be matched with (|4C)|1 to extract the n-graviton states. The (non-normalized) 
wave functional of the one-graviton state with momentum p and polarization e, for example, reads 



%, e [hij] = 6(h ) h^^p) ^ [h l3 ] 



(42) 



and so on. 



4 Newton potential from the propagation kernel 

In the presence of external static sources, the energy of the lowest energy state must be the Newton 
self-energy of the external source. Therefore the Newton self-energy can be extracted from the field 
propagation kernel as its lowest Fourier component in T. This procedure has proven effective in YM 
theory jlll . In particular, in the abelian case, the external source can be taken to be static, and 
the lowest energy state is characterized by the Coulomb self-energy. In the non-abelian case, on the 
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other hand, the external sources cannot be static, reflecting the exchange of colour charges between 
the external sources and the system. Inserting an external source in GR is a more delicate procedure, 
because the Newton potential emerges not only in the non-relativistic limit, as the Coulomb potential 
in YM theory, but also in the low gravity limit. This means that we can follow the same procedure one 
uses in YM theory, but we expect for the Newton potential to emerge only in the linear approximation. 
To introduce an external source S^v {%■> t) , we consider the lagrangian 



(43) 



The source is the densitized matter energy-momentum tensor J^ v = <?T M „. In the linear approx- 
imation around Minkowski space, the lagrangian P^l becomes 



(44) 



We characterize a static external source with the condition that only the component j7oo = P be 
different from zero. p(x, t) is thus the energy density of the source. Covariance is broken since the 
source itself defines a preferred frame. The source term in (|44J1 reads then p + h 00 p. Note that the 
conservation low V^T^y = at first order reads dop = 0, consistently with the fact that the source is 
static. In the temporal gauge ho^ = the Lagrangian becomes 



1 [hij , /iq m 



0; JfJi = C[hij] 



(45) 



Notice that we seem to lose the coupling between the gravitational field and the external source. This 
is analogous to what happens in YM, where the coupling term Aqp to an external static source is 
killed by the temporal gauge Aq = 0. This seems to prevent us from coupling external sources to the 
field in the temporal gauge. But the coupling can nonetheless be obtained, because part of the gauge 
degrees of freedom turn out to describe the source Below we show how this happens. Indeed, 
staring from lagrangian l|45|) and following the same steps as in the previous sections we arrive at the 
following expression for the propagation kernel, 



• exp ■ 



d 3 p 



2rp 2H ijk i(p) + H ijkl (p) 

P J ^2 L> %3 L) kl + 



H ijkl (p) H m (p) i p2 DijDiA 2pkpiD ^ 



+P 



2p 2 T 

Hjjki (p) cos pT - Hjju (p) 
2 sinpT 



DikDji+DuDjk-DuD 



ij^kl 



(46) 



where 



Hijkl {x, y) = h'lj {y)ti k \ (x) + h' t] (y)h' kl (x), 
Hijkl (2, V) = Kj [y) h' kl (x) + h\j (y) h' k \ (x) . 



(47) 



As before, the integrals over the e s implement the Wheeler-DeWitt constraint, but in presence of 
matter this is now didjh'fj — V 2 /i" + p, or, in momentum space, 2p 2 h — pip)- Using these constraints 
we get 



W[h>j,h»j,T] = M{T) exp { j J |^ 



smpT 



(48) 
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where 

Eo = ~ [ d?x [ dV-^ 1 ^-, ™ = I p{x) d*x. (49) 

32tt7 J \x-y\ J 

Therefore, the ground state is now weighted by the classical expression for the Newtonian self-energy 
of the external source and by the its rest mass. In the presence of sources, the field propagation nicely 
codes the Newtonian interaction. 



5 Summary 

We have studied the field propagation kernel of general relativity in the temporal gauge, W[g'ij, g'^T]. 
This is given in equation (|25[) . It can be interpreted as a 3-geometry to 3-geometry transition ampli- 
tude, or as a projector on the solutions of the quantum constraints. It is independent from T when g'^ 
and g'?j are defined on compact spaces; it depends on the asymptotic time if they are asymptotically 
flat. When g'^ and g' ^ ar e close to Sij, it can be computed explicitly to first order: the resulting 
expression is given in (|39|l . It is independent from the spin-one and spin-zero longitudinal compo- 
nents hi of the linearized field, and concentrated on the spin-zero component ho = values. From its 
form, the linearized vacuum, given in l|41|l . and n-graviton states, as in (|42|l . follow immediately. By 
coupling a static source to the field, we can extract the Newtonian self-energy (|49|l from the linearized 
expression of the kernel. 

These results are based on an accurate implementation of the FP procedure. The projection on 
the solutions of the constraints is implemented by the integration over the FP gauge parameters. This 
happens in a more subtle way than in YM theory, since temporal diffeomorphisms affect the boundary 
of the action. 

We have disregarded ultraviolet divergences. There are tentative background-independent def- 
initions of quantum gravity that define 3-geometry to 3-geometry transition amplitudes free from 
ultraviolet divergences. In order for these to have an acceptable low energy limit, they should reduce 
to the propagator kernel computed here, for values of their arguments close to flat space. 

Finally, showing that the background-independent definitions of quantum gravity lead to the 
correct Newtonian interaction at low energy has proven so far surprisingly elusive. Here we have 
shown how to extract this interaction from the field propagation kernel, by adding an external source. 



MT thanks the MIUR (Italian Minister of Instruction University and Research) and INFN (Italian 
National Institute for Nuclear Physics) for partial support. CR thanks the physics department of the 
University of Rome "La Sapienza" and FM, MT and SS thank the CPT for hospitality during the 
preparation of this work. 



Appendix A: Hamilton function of linearized GR 



In this appendix we evaluate the Hamilton function of linearized GR. This is the value of the grav- 
itational action on a classical solution with given boundary values, expressed as a function of these 
boundary values pQ. This will be used in the next appendix, to obtain the expression (|39|l of the 
propagation kernel. To simplify the calculations, we work with Euclidean signature, using analytic 
continuation to switch back to Lorentzian signature at the end of the calculations. In this appendix, 
x° is the Euclidean time, defined by x° = it. 

The quadratic part of the action of the linearized field, in the temporal gauge is 



d A x 



d k h ij d j h ik + d k h ij d i h jk 



d 0, h lJ d a h ij 



dih ij djh - d 3 h kj d k h + d a hd a h 



(50) 
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In the following calculations the linear term in (|31|1 plays no role. The equations of motion obtained 
by varying this action are 

□ h l3 - did m h m3 - d 3 d m h im + d l d j h + 5 io d m d n h mn ~ h = 0. (51) 

We are interested in evaluating the action on the classical solution with boundary data h'^ and h'^. 
Using the field equations and the boundary data, the action can be written as 



S cl = --(8*51- S ab 8 M ) I d A x 



h" ab (x)d*h kl (x) - ti ao {x)d%h kl {x) 



ah / 



x°=T 



x°=0 



(52) 



where the hij (x, t) is the classical solution interpolating between h'^ and h"j . This classical solutions 
can be constructed by means of a Green function Gij k i(x,y), defined by 

nG ijH {x,y) - d i d m G mjk i(x,y)-d j d m G imH (x,y) + d i d j G m mkl (x,y) + (53) 

+ S ijd ^G mnkl{ *,y) ~ S v BG- mkl (x,y) = \ {5lkh + S u S jk )S^x-y), 

with boundary conditions Gij k i(x,0;y) = Gij k i(x,T;y) = 0. To find it, we Fourier transform in the 
spacial components, 

G m (x,v)= J -^G ljkl {p-x\y°)e-^-^- 

then decompose Gij k i into covariant tensors and match the corresponding terms in (|53|l . We obtain 
GijM (p; x°, y°) = A(p 2 ,x°, y°) p l p J p k pi + 

+B(p 2 ,x°, y°) - piPjPkPi + p 2 (piPk8ji + PjPiS ik + PiPiS jk + PjP k Su - PiPjSki - PkpAj) 
+C(p 2 ,x°,y°) p A (S ik Sji + SuSjk - 8 i0 8 k i) + 

-p 2 {piPk8ji +PjPiS ik +PiPiS jk +Pjp k Su -piPj5 k i -p k pi8rj) +PiPjp k pi 



(54) 



where the functions A,B,C are defined as follows: 



B(p 2 ,x°,y°) = (y»-T 



12p 2 T 
2p 4 T 2p 4 



12p 2 



- (py° ~ pT 

2p 5 sinh pT ' 2p 5 sinh pT 



2 o Os sinh px° sinh (py° - pT) Q/ ^ i smh(px° - pT) sinhpy af _ Q 0n 

^\P jX ,y ) — — 



(x° - A 



The Green function allows us to write the solution of the field equations with given boundary data 



ti mn {x)d%G ijkl (x,y) 



c°=T 



(55) 



By inserting l|54|) into (|55|l and this into l|52|) . we obtain the Hamilton function of unconstrained 
linearized GR 



S c i[h' lv h'[ 3 ,T] - --J 



2 2H m {p) + H i]kl {p) 

— P-L 77; l->ijL>kl + 



Hjjki(p) - Hjjkijp) 
12 " ' " "' ' 2p 2 T 



p 2 DtjDki + PiPkDji + PjPiD lk + p l piD jk + Pjp k Du - 2piPjD k i - 2p k piDi 



+P 



Hjjki (p) coshpT - H ljkl (p) 
2 sinhpT 



I),, J), 1 ■ I) ,1 n,j) H 



(56) 
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where Dij and -Hy/u are given in 1(13(1 and (|47|l . 



Appendix B: Evaluation of the propagation kernel 

Here we illustrate the derivation of 1(39(1 from 1(37(1 . To perform the gaussian integral over hij , we write 
the field as the solution of the classical equations of motion h®j with the given boundary data, plus a 
fluctuation Qj . The boundary data are then 

&(0)=&(T) = 0. (57) 



Linear terms in Qj vanish because of the equations of motion, the integration over Qj yields a nor- 
malization factor J\f(T) depending of T, but independent of the boundary data. The remaining 
exponential, independent of the perturbation field, is the Hamilton function, computed above. Using 
this, we obtain for the Euclidean propagation kernel 



WEih'^h'lpT] = N(T) S {didjh'd ~ V 2 /i") S {didjti^ - V 2 /i 

1 f d 3 p 



J ^(t,x) ex P - \ J J# X J d s y e -m-*i 



P ^2 VijL) k i + 



p 2 D i3 D kl + 



2p 2 T 

+PiPkDji +PjPiD lk +PiPiD jk +pjp k D u - 2piPjD M - 2p k piDij 

. H tjki(x,y)coshpT - Ht jkl (x,y) , x 1 

+P — 2sinhpT [DikDji + D a D jk - D i:j D k ij S . (58) 

Here the apex e indicates a spatial diffeomorphism e l only on the field h'{j , 

Ht jkl {x,y) = h^l^U^ + h'^h'^x), (59) 
H! jkl (S,y) = h^WtiuW + h^m^Z). (60) 

The term in the last line of the expression 1)58(1 has the structure of the propagation kernel for a 
harmonic oscillator (see for instance [Q). 

The next step is to consider the integration over the spatial diffeomorphisms; the integrals can 
be performed separately for the longitudinal and the transversal parts 8 of e l . When we perform the 
integration over the transversal spatial diffeomorphisms, the result we obtain is that only the term 
with the structure of a harmonic oscillator survives. On the other hand, the integration over the 
longitudinal part c^e 1 gives again an implementation of the scalar constraint, in the form 

S (didjh'^ - V 2 h" - didjh^ + V 2 ti) . (61) 

The redundance of the integration over the e°'s and the die 1 can be understand as follows: in the 
hamiltonian formalism, the scalar constraint is proportional to the lapse function N. In the "Lapse=l, 



'We have the usual decomposition of a 3-vector, e l = + e T = ^jd l dje ] + e ! — -^ I d l djt 3 . 
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Shift=0" gauge we are working it, the transformation of the lapse under a diffeomorphism is given 
by (see for instance 9 ) N e — 1 + doe + die 1 . Therefore, both the integration over e° and over die 1 
provide variation of the lapse and consequentely an implementation of the scalar contraint. 

In the case of matter, the redundant implementation of the scalar constraint (|61|) . coming from 
the integration over the longitudinal spatial diffeomorphisms, is consistent with the other <5-functions, 
because the source is static and thus p{T) = p(0). 

The final expression can be extended to Lorentzian signature by means of the analytic continuation 
T — ► iT, and we obtain 



W\ti l3 , h% , T] = M(T) S (didjh'd ~ V 2 /i") 8 (didjh'ij - V 2 ti 
• exp 



£ / ' ^ f f d3p o-m-s) - ( f > ill cospT ~ & ill 

\lj J J (2^ e P 2smpT 

■ {D ik D jt + D u D jk - DtjDu) | , (62) 
which coincides with Ip^jl . once the spacial indices are contracted. 



Appendix C: Linearized Einstein equations in the temporal 
gauge 

For completeness, we write here the linearized Einstein equations in the temporal gauge ho^ = 0, 
extensively utilized in this paper. To do so it is convenient to work in the Fourier space with the 
Poincare-irreducible components introduced in (|32|l . These can be obtained from the field as follows 



ho 


— -jhkiDki, 


K = 


, PkPi 

- tiki 9 , 




= n-ki — Mi, 
P 


h TT 


= h M ID., 



ii - ^A^hJ • (63) 

We consider the case of a dust distribution with energy density p(x). So the only non-vanishing 
component of the energy-momentum tensor is Too = P- As pointed out in section Q, the source has 
to be static in the linearized case because of the continuity equation. We will found this result also in 
the study of the Einstein equations; this is not surprising because the continuity equation comes from 
the Bianchi identities, which are properties of the Einstein equations themselves. The vacuum case is 
easily reconstructed setting p = 0. 

The linearized equations for a dust distribution are 

S Q cP V " d»d a Ku ~ d v d a h afM + 8^h a a = | V (64) 

They can be written in the temporal gauge as follows 

h'i - ~ , (65) 
-d j h ij + d i h k k = 0, (66) 
nhij - did k h kj - dj&hji + didjh k k = fa , (67) 

where the dot indicates the derivative respect to the zero component and □ = d a d a . From the 
time-equations (|65I66|I we obtain the constraints in terms of the quantities defined in (|63(l 

h = 0, hf = 0, h\= 9 -. (68) 
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On the other hand, imposing the constraints on the motion equations l|67(l we have 

UhJ7 = 0, 2p 2 h = p. (69) 

The first is a wave equation for the traceless transverse components, which are the only physical 
degrees of freedom. The second is the Wheeler-deWitt equation in the Fourier space, as can be easily 
seen from the definition of hg in l|63l) . As it was previously said, the first of the constrains (|68|l together 
with the Wheeler-deWitt equation imposes the source to be static. 

In terms of the Poincare-irreducible field components, the quadratic term of the Lagrangian can 
be written as 

4£< 2 ' = hJJUhJJ + 2h Dh + Ahok - 2hJhJ. (70) 
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